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Abstract. By detailed first-principles calculations we show that the Fermi energy 
and the Rashba splitting in disordered ternary surface alloys Bi x Pbj / Sbi_ x _ 1/ /Ag(lll) 
can be independently tuned by choosing the concentrations x and y of Bi and Pb, 
respectively. The findings are explained by three fundamental mechanisms, namely 
the relaxation of the adatoms, the strength of the atomic spin-orbit coupling, and 
q ' band filling. By mapping the Rashba characteristics, i.e. the splitting &r and the 

O . Rashba energy £"r, and the Fermi energy of the surface states in the complete range of 

concentrations, we find that these quantities depend monotonically on x and y, with 
a very few exceptions. Our results suggest to investigate experimentally effects which 

J^ , rely on the Rashba spin-orbit coupling in dependence on spin-orbit splitting and band 

\Q '. filling. 
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1. Introduction 

In the emerging field of spin electronics, proposed device applications often utilize the 
Rashba effect jl] in a two-dimensional electron gas (2DEG). A prominent example is the 
spin field-effect transistor [2j in which the spin-orbit (SO) interaction in the 2DEG is 
controlled via a gate voltage [HIS]- Other examples are a high critical superconducting 
temperature which shows up in materials with a sizable spin-orbit interaction [5] and 
the spin Hall effect [BHH]. 

The Rashba effect relies on breaking the inversion symmetry of the system and, 
consequently, shows up in semiconductor heterostructures and at surfaces. The breaking 
of the inversion symmetry results — via the spin-orbit coupling — in a splitting in the 
dispersion relation of electronic states which are confined to the interface [lj. In a 
simple model for a two-dimensional electron gas, a potential in z direction confines 
the electrons to the xy plane. The Hamiltonian of the spin-orbit coupling can thus be 
written as 

H B0 = jR.(cr x dy - a y d x ), (1) 

where the strength of the SO interaction is quantified by the Rashba parameter tr. 
Employing a plane-wave ansatz yields the dispersion relation 

h 2 k 2 
£ ± (£|])=£oH-^± 7R |% (2) 

where m* is the effective electron mass. The split electronic states are labeled by + 
and — ; their spins lie within the xy plane, are aligned in opposite directions, and are 
perpendicular to the wave vector k\\. 

In a real system, the Rashba parameter 7r comprises effectively two contributions 
[TO] . The 'atomic' contribution is due to the strong potential of the ions (atomic 
spin-orbit coupling). The 'confinement' contribution is due to the structural inversion 
asymmetry which can be viewed as the gradient of the confinement potential in z 
direction. The larger this gradient and the atomic spin-orbit parameter, the larger 
7r and the splitting 

, |ra*|7 R 



m 



(3) 



which is defined as the shift of the band extremum off the Brillouin zone center [k\\ = 0). 
Another quantification of the splitting is the Rashba energy 

h 2 kl m*7 R 

Er = ~ 2^ = " 2/iT' (4) 

that is the energy of the band extremum with respect to the energy E for which the 
bands cross at k\\ = 0. 

The above dispersion relation suggests to distinguish two energy ranges. Region I 
is defined as the energy range between E Q and the band extrema (E G [E — E-&, E ] 
for positive m* or E G [E ,E + E-&] for negative m*) [11]. Region II comprises the 
other range of band energies (E > Eq for positive m* or E < Eq for negative m*). The 
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density of states in region I is singular at the band extrema and decreases towards E 
while in region II it is constant. 

In view of designing device applications and investigating fundamental effects, it is 
desirable to tune both the strength 7^ of the Rashba spin-orbit coupling and the Fermi 
energy Ep of the 2DEG. In a semiconductor heterostructure, this can be achieved by an 
external gate voltage and by doping of the semiconductor host materials. At a surface, 
these quantities can be affected by adsorption of adatoms [T2|fT3] , by surface alloying jTT] , 
and by changing the thickness of buffer layers (e. g. in Bi/(Ag) n /Si(lll) [14J). Recently, 
a ferroelectric control has been proposed |15j . 

Surface states in surface alloys show an unmatched Rashba splitting [IE], as 
has been investigated in detail by scanning tunneling microscopy as well as by spin- 
and angle-resolved photoelectron spectroscopy. They are convenient systems for 
testing fundamental Rashba-based effects. The ordered surface alloys Bi/Ag(lll), 
Pb/Ag(lll), and Sb/Ag(lll) have been investigated by first-principles calculations 
and in experiments [IEHT8]. These three systems differ with respect to their Rashba 
characteristics k R and E R , and by E . The challenge we are dealing with is how to tune 
these properties independently. 

The basic idea is as follows. Bi/Ag(lll) has a large splitting and occupied sp z 
surface states, while Pb/Ag(lll) has a large splitting and unoccupied sp z surface states. 
In a disordered binary alloy Bi x Pbi_ x /Ag(lll) the Fermi energy can be tuned by the 
concentration x, while keeping a large spin splitting. In contrast, Sb/Ag(lll) has 
occupied surface states with almost the same binding energy as those in Bi/Ag(lll) 
but a minor splitting. This allows to tune mainly the spin splitting but keeping the 
Fermi energy in Bi x Sbi_ z /Ag(lll). Thus, by an appropriate choice of concentrations 
x and y in a ternary alloy Bi x Pb ?/ Sbi_ x ._ ?/ /Ag(lll) we expect to tune the Fermi energy 
and the splitting independently. In particular, one could access the region I between E 
and the band maxima which is important for high-temperature superconductivity [5]. 

We report on a first-principles investigation of disordered surface alloys 
Bi x Pb ?/ Sbi_ :r _ l/ /Ag(lll), performed along the successful line of our previous works on 
both ordered and disordered alloys [TT],[T6|,[T8] . Since all ordered and disordered binary 
alloys show a v^ x v^3-R30° surface reconstruction, we assume this geometry also for the 
ternary alloys. The resulting substitutional disorder is described within the coherent 
potential approximation. 

The paper is organized as follows. Our computational approach is sketched in 
section [2j The results are discussed in section [31 for binary alloys in section 13.41 and for 
ternary alloys in section 13.51 We give conclusions in section HI 

2. Computational aspects 

We rely on our successful multi-code approach, based on the local density approximation 
to density functional theory. Because this is described in detail elsewhere [15], we 
deliberately sketch it in this paper. 
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The surface relaxations of ordered surface alloys were determined using the Vienna 
Ab-initio Simulation Package (VASP) |19j . well-known for providing precise total 
energies and forces. The relaxed structural parameters serve as input for first-principles 
multiple-scattering calculations. Our Korringa-Kohn-Rostoker (KKR) method already 
proved successful for relativistic electronic-structure computations of Rashba systems 

[11120]. 

The central quantity in multiple-scattering theory is the Green function [21] 

G(f n , C; E,k)=J2 Z A(r n , E)r^(E, k)Z^ m , Ef 

AA' 

-S™Y,m?<\E)JW>\E)\ (5) 

A 

where Z and J are regular and irregular scattering solutions of sites n and m at energy 
E and wavevector k, respectively. f n is taken with respect to the position R n of site 
n{f n = f— R n ). r < (r>) is the lesser (larger) of r n and r' n . A = («,//) comprises 
the relativistic spin- angular-momentum quantum numbers [21J. The scattering-path 
operator r is obtained in standard KKR from the so-called KKR equation [21], or in 
layer-KKR from the Dyson equation for the Green function [22J. 

The local electronic structure is analyzed in terms of the spectral density 

N n (E;k) = --3Tr G(f n , f n , E, k). (6) 

n 

By taking appropriate decompositions of the trace, the spectral density provides 
information on spin polarization and orbital composition of the electronic states. 

Substitutional ternary alloys BL r Pbj / Sbi_, r _ 2/ /Ag(lll) are described within the 
coherent potential approximation (KKR-CPA), in which short-range order is neglected. 
From the agreement of the theoretical data with their experimental counterparts for 
the binary alloys Bi x Pbi_ z /Ag(lll) [TT], we conclude that short-range order is of minor 
importance in these systems. Hence, we applied the KKR-CPA also for the ternary 
alloys. 

The effect of the disorder can be understood as a self-energy [23]. As a consequence, 
the spectral density of the disordered alloys becomes blurred (or smeared out) as 
compared to that of the ordered alloys. 

3. Results and discussion 

3.1. Geometry 

Relaxations have been determined by VASP for the ordered alloys, with v3 x v3i?30° 
reconstruction and face-centered-cubic (fee) stacking (VASP cannot treat substitutional 
disorder within the CPA). It turns out that the relaxations of Sb, Bi, and Pb are in 
accord with their atomic radii. To be more precise, the outward relaxations are 9.6, 
15, and 18 percent of the Ag(lll) bulk interlayer spacing (2.33 A), respectively, with 
respect to the positions of the Ag atoms in the topmost layer. Being negligibly small, 
in-plane displacements of Ag atoms are not considered. 
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Since all ordered and disordered binary alloys show a \/3 x ^^30° surface 
reconstruction, we assume this geometry also for the ternary alloys. The relaxations 
of the disordered surface alloys were linearly interpolated, in dependence on the 
concentrations of the constituting elements Bi, Pb, and Sb. This assumption is within 
the spirit of the CPA; being a mean-field theory, a disordered system is described by an 
effective medium. Likewise the relaxation should be taken as a concentration-weighted 
average. We are aware, however, that in real samples, the relaxations of the constituting 
individual atoms could differ, as might be checked by scanning tunneling microscopy. 

3.2. Mechanisms which influence the Rashba-split surface states 

Before presenting details of our calculations, a brief discussion of the general trends and 
mechanisms is in order. For tuning the Fermi energy and spin splitting independently, 
the underlying mechanisms should be independent as well. 

A first mechanism is relaxation. The outward relaxations of Sb, Pb, and Bi are 
in accord with their atomic radii; the larger the atomic radius, the larger the outward 
relaxation. The relaxation is accompanied by a charge transfer from the atomic sphere 
to the surrounding: the larger the relaxation, the larger the charge transfer [21]. This 
mechanism determines the energy position of the degenerate point Eq — cf. (J2]) — and, 
consequently, the Fermi energy or band-filling of the surface states (2DEG). 

A second mechanism is the atomic spin-orbit parameter. Bi and Pb are heavy 
elements with large SO parameter (1.25 eV for Bi and 0.91 eV for Pb [25]), in contrast 
to the lighter element Sb (0.4 eV [25]). The Rashba splitting depends both on the atomic 
SO-coupling strength and the potential gradient [TD]. Since the latter should not differ 
considerably among the considered systems, the spin splitting is mainly determined by 
the atomic SO coupling. We expect that with increasing Sb content, the spin splitting 
decreases. 

A third mechanism is electron doping or band filling. Pb has one electron less than 
Bi (Zpb = 82, Zb[ = 83). Within a rigid-band model, the surface states in Pb/Ag(lll) 
are shifted to higher energies, as compared to those in Bi/Ag(lll). This picture is 
confirmed by experiments and first-principles calculations |llj . 

3.3. Ordered surface alloys 

The ordered surface alloys Bi/Ag(lll), Pb/Ag(lll), and Sb/Ag(lll) have been studied 
previously in detail [T6HT8"] . They show two sets of surface states; a first set is unoccupied 
and consists mainly of p x p y orbitals (for Bi/Cu(lll), see [2S]). In this paper, we focus 
on the other set which is either completely or partially occupied and consists of sp z 
orbitals. The effective mass m* of both sets is negative, implying a negative dispersion. 

Sb/Ag(lll). We address briefly the abovementioned relaxation mechanism by 
considering two cases for Sb/Ag(lll): (i) an Sb relaxation as calculated by VASP 
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(9.6 %) and (ii) an artificial relaxation of 25 %. The charge transfer from the Sb muffin- 
tin spheres to the surrounding is increased for the larger relaxation (2.05 % as compared 
to 0.94 %, with respect to the nominal valence charge; cf. [21]). Consequently, the surface 
states are shifted towards higher energies by 0.16 eV, as obtained from the degeneracy 
point Eq. Further, the spin splitting &r becomes increased as well (0.03 A -1 as compared 
to 0.02 A -1 ). This corroborates that the relaxation mainly affects the crossing point E 
(or Fermi energy) rather than the spin splitting. 

3.4- Disordered binary alloys 

Bi x Pb\^ x /Ag(l 1 1) . In the disordered binary alloy Bi x Pbi_a;/Ag(lll), which has been 
studied previously [TT], the ratio of the Rashba energy Er and the Fermi energy Ep can 
be chosen within a wide range, in dependence on the Bi concentration x. For both Bi 
and Pb, 0.99 % of the atomic charge atom is removed from the muffin-tin sphere, which 
is in agreement with the close outward relaxation of Bi and Pb (15% and 18%). As 
noted before, Pb has one valence electron less than Bi, which explains the sizable shift 
of the surface states to higher energies (band-filling mechanism; cf. the panels on the 
right-hand side of figured]). Although the relaxation is of the same order, the splitting 
is smaller for Pb (topmost panel in figure [T|). This can be attributed to the smaller 
atomic spin-orbit parameter of Pb (0.91 eV for Pb and 1.25 eV for Bi [25]). 

Bi x Sbi_ x /Ag(lll). Recently, the surface states of the disordered binary alloys 
Bi x Sbi_ :r /Ag(lll) were mapped out by angle-resolved photoelectron spectroscopy. The 
momentum offset &r evolves continuously with increasing Bi concentration x. The 
splitting decreases sizably for x < 0.50 [27]. 

In theory, the outward relaxation of Bi is larger than for Sb (15 % and 9.6 %, 
respectively). Consequently the charge which is removed from the Sb sphere (0.94%) is 
smaller than that of Bi (0.99%). Since Bi and Sb are iso-electronic, with valence-shell 
configuration 5p 3 and 6p 3 , Eq remains almost unaffected by x, as can be seen in the 
bottom row of figure [1] The spin splitting for Sb is much less than for Bi, in agreement 
with the atomic spin-orbit parameter (0.4 eV and 1.25 eV). In accord with experimental 
results, the Rashba splitting k^ evolves with Bi concentration x. 

To elucidate further the effect of the relaxation, we calculated the splitting 
of Bi .6Sbo.4/Ag(lll) for two relaxations. The interpolated relaxation for 

Bi . 6 Sbo.4/Ag(lll) is 12.8% (shown at (x,y,z) = (0.6,0.4,0.0) in figured]), for the 
artificial relaxed system the outward relaxation is taken as 19 % (not shown here). The 
charge transfer for the two systems is very close, and the difference in the splitting is 
negligibly small. Hence, the splitting is negligibly sensitive to the relaxation, as was 
already established for Sb/Ag(lll). 

PbySbi_y/Ag(lll). To complete the picture of the binary alloys we turn to 
Pb y Sbi_ y /Ag(lll), for which experimental results are not available. The trends which 
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Figure 1. Surface states of disordered ternary alloys Bi x Pb 2/ Sbi_ a; _j / /Ag(lll) along 
T-K of the two-dimensional Brillouin zone. The spectral density at a heavy-element 
site Hi x Pb y Sbi- x -y is depicted as linear gray scale, with dark gray corresponding to 
high spectral weight; cf. ((6]). 



Ternary surface alloys 8 

have been discussed before are as well found in these alloys (cf. the panels on the 
left-hand side of figure [1]). As the Pb concentration increases, E shifts down from 
Ep + 0.6 eV to Ep — 0.4 eV, implying that the surface states become completely filled 
at about y = 0.3. As for Bi x Sbi_ x /Ag(lll), the spin splitting increases with y. 

3.5. Disordered ternary alloys Bi x Pb y Sbi_ x _ y /Ag(lll) 

Having established the ingredients which are necessary for independently tuning 
the Fermi energy and the spin splitting in the surface alloys — by investigating the 
disordered binary surface alloys — we now mix them to disordered ternary alloys 
Bi x Pbj / Sbi_ :r _ l/ /Ag(lll). By choosing appropriate concentrations x and y, the 
degeneracy point E and the Rashba splitting are tuned. Note that the splitting /cr 
and the Rashba energy E-& are not fully independent; both can be expressed (in a free- 
electron model) in terms of the effective electron mass and the Rashba parameter [cf. 
(ED and gl)]. 

In figured] the surface-state dispersions of ternary alloys BL c Pbj / Sbi_. r _ y /Ag(lll) 
are shown. The concentrations x and y have been varied in steps of 0.2. A common 
feature of the spectral density of the binary and ternary alloys is a finite lifetime of the 
spectral density, which is the consequence of the substitutional disorder. 

The Rashba characteristic of the ternary alloys follow the general trends of the 
binary alloys which have been discussed before. In the ternary alloys with larger outward 
relaxation (i.e. the Bi- and Pb-rich compounds), the degenerate point Eq shifts toward 
higher energies (main mechanism: relaxation). The larger the concentration of heavy 
elements Bi and Pb as compared to the Sb concentration, the larger the splitting k^ 
(main mechanism: atomic spin-orbit parameter). The degenerate point E shifts upward 
with increasing Pb concentration (main mechanism: band filling). 

The shift k^ of the surface states in reciprocal space versus concentrations x and y 
is shown in figure [2] (top). As expected, the smallest splitting (dark blue) shows up for 
Sb/Ag(lll) (z = 1 — x — y = 1), while the largest (dark red) corresponds to Bi/Ag(lll) 
(x = 1). For Pb/Ag(lll), & R is of intermediate order (green/yellow). Surprisingly, the 
splitting is not monotonic, as one might have expected in a rigid-band picture. For 
example, &r shows a local minimum at (x, y, z) w (0.4, 0.4, 0.2). 

As &r, the Rashba energy E^ depends monotonously in a large range of 
concentrations (bottom in figure [2]). Sizable Rashba energies are found mainly for Bi- 
rich alloys, say for x > 0.5. This implies that for accessing region I, Bi-rich surface 
alloys are inevitable. For smaller x (blue areas in the bottom panel of figure [2]), the 
energy range of region I could be too small to be employed in experiments. 

The energy E of the degeneracy point depends almost linearly on the heavy 
elements' concentrations x and y (figure [3]) . For equal Bi and Sb concentrations (x = z) 
it is nearly constant; upon adding Pb, E shifts up. For systems with about 40 % of Pb 
concentration, Eq is very close to the Fermi level Ep, so that the latter lies in region 

i m 
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In summary, the above results support that both Fermi energy and spin splitting in 
the surface states can be tuned independently, as is readily apparent from the different 
shapes in figure |2] and figure [3j A very interesting region in the ternary plots is around 
(x,y,z) ~ (0.6,0.3,0.1), where the degenerate point E and the Fermi energy E F 
coincide. Keeping the Sb concentration constant and changing the Pb concentration 
of about 10 % is accompanied by transition between region I and region II, while &r 
and £r are almost constant. It is also possible to tune Er and /cr while keeping the 
position of degenerate point constant. The changes of &r and £r are not independent 
but /cr depends more sensitive on the concentrations than Er. 

4. Conclusions 

Disordered ternary surface alloys Bi a; Pb 3/ Sbi_ a ;_j,/Ag(lll) allow to fabricate a two- 
dimensional electron gas with specific Rashba spin-orbit splitting and Fermi energy 




Figure 2. Spin splitting in disordered ternary alloys Bi x Pbj,Sbi_ 2; _j,/Ag(lll). Top: 
The surface-state displacement /cr (in reciprocal space) is depicted as color scale as a 
function of Bi concentration x, Pb concentration y, and Sb concentration z = l—x — y. 
The color bar on the right is in units of A -1 . Bottom: Same as in the top but for the 
Rashba energy E-r. The color bar is in eV. 
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Figure 3. Surface-state energy in disordered ternary alloys Bi a; Pbj,Sbi_ K - a /Ag(lll). 
The degeneracy energy E of the surface state, with respect to the Fermi level E F , is 
depicted as color scale as a function of Bi concentration x, Pb concentration y, and Sb 
concentration z = 1 — x — y. The color bar on the right is in eV. At negative energies, 
the surface states are fully occupied (blue area). 



which can be investigated by surface- scientific methods (scanning tunneling probes and 
especially photoelectron spectroscopy). In particular, the important transition from 
energy region I (that is, the Fermi energy Ep lies above the degeneracy point E ) to 
region II (E F below E ) can be studied for different strengths of the Rashba spin-orbit 
coupling. Thus, the present study may stimulate further experiments on Rashba systems 
and their unique properties. 
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